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ABSTRACT 

Context. The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission 
measure. The cross-sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the 
heating rate. 

Aims. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation 
in terms of steady heating models. 

Methods. We study the area response with a time-dependent 2D magnetohydrodynamic (MHD) loop model, including the description 
of the expanding magnetic held, coronal heating and losses by thermal conduction, and radiation from optically thin plasma. We run 
a simulation for a loop 50 Mm long and quasi- statically heated to about 4 MK. 

Results. We find that the area can change substantially with the quasi-steady heating rate, e.g., by ~40% at 0.5 MK as the loop tem- 
perature varies between 1 MK and 4 MK, and, therefore, affects the interpretation of the differential emission measure vs. temperature 
(DEM(T)) curves. 
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1. Introduction 

An important challenge in trying to understand the structure 
and heating of solar active regions is to explain the relative 
amounts of material at different temperatures. This is typically 
expressed in terms of the differential emission measure distribu- 
tion (DEM(T)). Models that use steady heating and a flux tube 
geometry that slowly expands with height are unable to repro- 
duce the observed distributions. However, realistic flux tubes are 
expected to expand very rapidly with height between the chro- 
mosphere and corona, where the plasma beta (J3 = 8 nP/B 2 ) 
changes from being much greater than to much less than unity. 
The thin transition region where temperatures increase from 
chromospheric to coronal values occurs somewhere within this 
rapidly expanding throat. Thus, we expect the cross-sectional 
area of flux tubes to be a strong function of temperature. If the 
dependence has the right form, steady heating models may in 
fact be able to explain the observed DEM(T). This was demon- 
strated by Warren et al. (2010b) for the core of an active region 
(the generally hot central part of the active region that connects 
opposite polarity moss regions, which are the transition region 
footpoints of the hot loops). 

A potential difficulty with this picture has been pointed out 
by Tripathi et al. (2010). The height of the transition region 
depends on the coronal pressure and so moves up and down 
in response to changes in the coronal heating rate, either fast 
or slow. If the heating rate increases, so too does the pressure, 
and the transition region is forced downward. If the heating rate 

* The movie associated to Fig. 4 is available in electronic form at 
http : //www . aanda . org 


decreases, some of the overlying pressure is relieved, which al- 
lows the transition region to move upward. For gradual heating 
variations, the amount of vertical displacement is given approx- 
imately by 

where H g is the gravitational scale height of the chromosphere 
(approximately 500 km), g is the new heating rate, and go is the 
original heating rate. We have made use of the fact that g oc p 1 /6 
for quasi-static evolution, where p is the pressure. See Klimchuk 
(2006) for a similar expression appropriate to impulsive heating. 
A heating rate increase of a factor of 100, which increases the 
loop apex temperature by a factor of 3.7, causes the transition 
region to be displaced downward by approximately 2700 km. 
This is not small compared to the length scale of expansion 
in the flux tube throat (Gabriel 1976; Athay 1981). Thus, the 
cross-sectional area at any given temperature A(T) could change 
dramatically. This would alter DEM(T) so that a model is no 
longer compatible with observations. The very sensitive depen- 
dence of the A(T) relationship on the heating rate suggests that 
too much fine tuning is required for a steady heating model to 
have widespread applicability. 

The last statement involves a big caveat, however. It assumes 
that flux tubes remain rigid as the heating and pressure change. 
This is not the case in a real magnetohydrodynamic (MHD) sys- 
tem. With increasing pressure, a flux tube will expand laterally 
at the same time that the transition region is pushed downward. 
Hence, the area at a given temperature might remain constant, 
whereas it would decrease in a rigid tube. It is possible that the 
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A(T) relationship in the transition region could remain largely 
unchanged. (We note, however, that any lateral expansion would 
be even more pronounced in the corona, so the emission mea- 
sure relationship between the transition region and corona would 
then change.) The purpose of our study is to determine how A(T) 
varies in response to changes in the coronal heating rate, espe- 
cially slow changes corresponding to quasi-static evolution. If 
there is a strong dependence on the heating rate, then the steady 
heating interpretation of DEM(T) should be reexamined. 

In this paper we present a 2D MHD loop model that natu- 
rally accounts for loop expansion through the transition region. 
We have used this model to study the plasma response to the 
heating into and around the moss regions. At variance with ID 
models having a specified expanding cross section (Emslie et al. 
1992; Mikic et al. 2013), we use a full MHD description, i.e., 
the full set of MHD equations is solved together in a 2D spa- 
tial domain; this allows us to describe a beta-changing system 
and to have feedback between plasma and magnetic field in the 
critical region, i.e., around the transition region. As a conse- 
quence, and at variance with previous works, in our model the 
loop cross-sectional area is also a function of time and changes 
as the heating changes. Our 2D cylindrical coordinates make the 
description realistic for a 3D structure with cylindrical symme- 
try, which is a good approximation of a coronal loop (but see 
Malanushenko & Schrijver 2013 for deviations from circular 
loop cross sections). We have focussed on the expansion around 
the transition region, not in the whole loop (Peter & Bingert 
2012). Our numerical experiment is devoted to answering a very 
specific question on the interaction between the coronal heating 
and the cross section of the transition region. It is complemen- 
tary to self-consistent 3D MHD descriptions of coronal boxes 
(Bingert & Peter 2013). 

In Sect. 2 we describe our model, the related equations and 
numerical code, the initial and boundary conditions, and how we 
set up the simulation of interest. The simulation and its results 
are illustrated in Sect. 3 and discussed in Sect. 4. 


2. The model 

The concept here is to study the reaction and readjustments of a 
magnetic loop and of the plasma confined therein under the ef- 
fect of a gradual heating release inside it. We focus on a single 
loop and its surroundings, and therefore our approach is com- 
plementary to large-scale modeling (e.g., Carlsson et al. 2010; 
Gudiksen et al. 2011; Bingert & Peter 2011; Martinez- Sy kora 
et al. 201 la, b). 

We describe a coronal loop as a magnetic flux tube linking 
two far locations of the solar chromosphere. The footpoints are 
so far from each other that we can assume that they are rooted 
in independent chromospheres. For simplicity of modeling, we 
then describe the loop as straightened into a vertical flux tube 
linking two chromospheric layers at opposite extremes of the ge- 
ometric domain (top and bottom boundaries). We keep the grav- 
ity proper of a curved flux tube. The geometry of our domain 
is 2D cylindrical (r - z). The spatial domain is much broader 
than the cross section of a typical loop. The proper loop corona 
forms as soon as we put a heating excess in the region around 
the central axis (r = 0) and not elsewhere. The heating is trans- 
ported along the magnetic field lines and makes plasma expand 
from the chromospheres upwards, filling the space between the 
footpoints in the central region. 

Our model considers the time-dependent MHD equations 
including gravity (for a curved loop), thermal conduction 


(including the effects of heat flux saturation), a heating function, 
and radiative losses from optically thin plasma. 

The MHD equations are solved in the non-dimensional con- 
servative form, 
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are the total pressure and total energy density (internal energy 6, 
kinetic energy, and magnetic energy), respectively; t is the time; 
p = pmunu is the mass density; p = 1.265 is the mean atomic 
mass (assuming metal abundance of solar values; Anders & 
Grevesse 1989); mu is the mass of hydrogen atom; nn is the 
hydrogen number density; u is the plasma velocity; g is the 
gravity acceleration vector for a curved loop; I is the iden- 
tity tensor; T is the temperature; F c is the thermal conductive 
flux (see Eqs. (9)-(12)), the subscripts || and _L denote, respec- 
tively, the parallel and normal components to the magnetic field; 
k\\ = K\\T 512 andk ± = K ± p 2 /(B 2 T l/2 ) are the thermal conduction 
coefficients along and across the field; K\\ and K ± are constants; 
Ci S o is the isothermal sound speed; 0 = 1 is a free parameter and 
E sat is the maximum flux magnitude in the direction of 7^; A ( T ) 
represents the optically thin radiative losses per unit emission 
measure derived from the CHIANTI v. 7.0 database (Dere et al. 
1997; Young et al. 2003; Reale & Landi 2012) assuming coro- 
nal metal abundances (Feldman 1992) and H (s , t) is a function 
of space and time describing the phenomenological heating rate 
(see Sect. 3). We use the ideal gas law, P = (y - 1 )pe. We assume 
negligible viscosity and resistivity, except for those intrinsic in 
the numerical scheme. 

The calculations are performed using the PLUTO code 
(Mignone et al. 2007, 2012), a modular, Godunov-type code for 
astrophysical plasmas. The code provides a multiphysics, algo- 
rithmic modular environment particularly oriented toward the 
treatment of astrophysical flows in the presence of discontinu- 
ities as in the case treated here. The code is designed to make 
efficient use of massive parallel computers using the message- 
passing interface (MPI) library for interprocessor communica- 
tions. The MHD equations are solved using the MHD module 
available in PLUTO, configured to compute intercell fluxes with 
the Harten-Lax- Van Leer approximate Riemann solver, while the 
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second order in time is achieved using a Runge-Kutta scheme. A 
Van Leer limiter for the primitive variables is used. The evolu- 
tion of the magnetic field is carried out adopting the constrained 
transport approach (Balsara & Spicer 1999) that maintains the 
solenoidal condition (V • B = 0) at machine accuracy. 

The PLUTO code includes optically thin radiative losses in 
a fractional step formalism (Mignone et al. 2007), which pre- 
serves the second-order time accuracy, since the advection and 
source steps are at least second-order accurate; the radiative 
loss A values are computed at the temperature of interest us- 
ing a table lookup/interpolation method. The thermal conduc- 
tion is treated separately from advection terms through opera- 
tor splitting. In particular, we adopted the super- time- stepping 
technique (Alexiades et al. 1996) which has been proved to be 
very effective at speeding up explicit time- stepping schemes for 
parabolic problems. This approach is crucial when high values 
of plasma temperature are reached (e.g., during flares), explicit 
scheme being subject to a rather restrictive stability condition 
(i.e., At < (Ax) 2 /2rj, where r/ is the maximum diffusion coeffi- 
cient), because the thermal conduction timescale r con d is shorter 
than the dynamical one Td yn (e.g., Orlando et al. 2005, 2008). 

2.1. The loop setup 

We address a typical active region loop, with half length L = 
3 x 10 9 cm and temperature T - 3 x 10 6 K. 

We choose as the initial plasma condition a plane-parallel 
loop atmosphere with the temperature and density profile from 
a model of a hydrostatic loop (Serio et al. 1981) with an apex 
temperature about 8 x 10 5 K. For the chromospheric part of the 
loop we use a hydrostatic atmosphere with a uniform tempera- 
ture (i.e., 10 4 K). The temperature spans from 10 4 K in the chro- 
mosphere to 8 x 10 5 K in the corona. The density spans from 
-10 14 cm -3 in the chromosphere to -10 8 cm -3 in the corona. 
This atmosphere is much more tenuous and cool than the subse- 
quently heated one. The corona and the chromosphere are con- 
nected through a thin transition region where the temperature 
jumps from 10 4 K to 10 6 K in less than 10 8 cm. 

The computational domain (see Fig. 1) extends from -3.1 x 
10 9 cm to 3.1 x 10 9 cm in the z-direction (i.e., along the loop 
axis), and from r = 0 to r = 3.5 x 10 9 cm in the r-direction 
(i.e., across the loop). High resolution (dr - dz ~ 3 x 10 6 cm) is 
needed to describe appropriately the thin transition region from 
the chromosphere to the corona. The grid that we use is non- 
uniform but fixed. Figure 1 shows the pixel map in the domain. 
We have the maximum resolution around the transition region 
(\z\ « 2.4 x 10 9 cm, i.e., the layers where the colors change 
abruptly in Fig. 1) in the vertical direction and along the loop 
axis in the r- direction (r « 0). 

A fundamental ingredient of our study is that the loop where 
we inject the heating expands from the chromosphere to the 
corona. We want the magnetic field to have this expansion right 
from the beginning of the heating release. Therefore, our strategy 
involves two steps. First, we build the magnetic configuration in 
a tenuous, static and plane parallel atmosphere, i.e., we build an 
empty loop. Second, we inject the heating in a localized region 
of the domain to switch the loop on. For the first step, we start 
with a magnetic field that is in the vertical direction in the entire 
spatial domain, i.e., the field lines are straight and perpendicu- 
lar to the two dense independent chromospheres. The magnetic 
field and the plasma density are more intense around the central 
axis, i.e., around r = 0. This configuration is not at equilibrium, 
and a preliminary simulation is devoted to letting the configu- 
ration relax to equilibrium. The central initial density is chosen 



- 3 - 2-10 1 2 3 



r [10 9 cm] 

Fig. 1 . Top: color maps in logarithmic scale of density (left, [cm -3 ]) and 
temperature (right, [K]) before the loop heating is switched on. The 
chromosphere is red in the density map and blue in the temperature 
map. The magnetic field lines are also marked (white lines). Bottom: 
pixel map in the computational domain. Each box encloses 10x10 grid 
points. 


so that at equilibrium it is approximately equal to that of the 
surroundings. 

In the tenuous corona, around the central axis the magnetic 
field pressure is much higher than the plasma pressure, i.e., the 
ratio of thermal pressure and magnetic pressure is /5 <3C 1. As 
a result, the magnetic field expands considerably in the corona. 
The magnetic pressure is instead only a small perturbation in 
the chromosphere (J3 » 1) and is left almost unchanged there. 
We let the whole system evolve freely, keeping a small heating 
that compensates for radiative and conductive losses. A quasi- 
steady state is reached after t - 4000 s (the Alfven and sound 
crossing time are -50 s and -1000 s, respectively). The result 
of this preliminary step is a steady and stable loop atmosphere, 
with the strong expansion of the magnetic field lines from the 
chromosphere to the corona around the loop central axis, as we 
supposed (Fig. 1). The magnetic field intensity decreases from 
a few hundred G in the chromosphere to a few G in the corona 
(sufficient to confine the loop plasma). As shown in Fig. 1, the 
area expansion is a factor of 6 from the bottom of the chromo- 
sphere to the top of the transition region, another factor of 2 in 
the first 3000 km above the transition region, and another factor 
of 2 to the top of the loop (middle of the domain). 

The plasma atmosphere also readjusts to a new equilibrium: 
in Fig. 1 we see that in the magnetically-expanded configuration 
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the atmosphere is no longer strictly plane parallel both in the 
temperature and in the density, but the inhomogeneities are small 
and leave the corona very tenuous and much cooler than 1 MK. 
We also see some ripples in the chromosphere that do not affect 
the evolution when the loop heating is on. All plasma motions 
are a few km/s at most. 

Many different initial conditions have been tested, finding 
that the resolution in the transition region is critical for accuracy. 
We adopt axisymmetric boundary conditions along the symme- 
try axis (i.e., at r = 0), reflective conditions at the outer bound- 
ary (r = 3.5 x 10 9 cm) and reflective conditions, but with reverse 
sign, for the magnetic field components at the lower and upper 
boundary (z = ±3.1 x 10 9 cm). 


3. The loop simulation 

We use the results of the preliminary simulation as the starting 
point of a new simulation that studies the response of the loop 
to heating. Our scope is to study what level of variability is ex- 
pected in the moss at the base of a high-pressure loop. We check, 
for instance, whether the loop heating drives variations in the 
loop and in particular of its cross section in the transition re- 
gion, by considering the situation in which the loop is heated as 
smoothly as possible (e.g., Winebarger et al. 201 1; Warren et al. 
2010a, b). 

The heating function that we use is divided into two terms: 
the first is the static term H s that covers the entire corona uni- 
formly, and the second H t that slowly varies with time and cov- 
ers only a part of the loop: 


H(r,z,t) = H s + H t (r,z,t). 


(13) 


According to the scaling law (Rosner et al. 1978), the static term 
is set to H s = 4.2 x 10 -5 erg cm -3 s _1 , which is sufficient to 
sustain a static loop with an apex temperature of about 8 x 10 5 K 
(see Sect. 2.1). 

The second term varies with time and space as 


H t (s,t) = H 0 Mr)f z (z)Mt), 
where 

fr(r) = e -r2/2cr », 

( 0 z < -zo 

fz(z) = \ 1 -Zo<Z<Zo 

[ 0 z > zo, 

f(t) = i° t<to 

JA ’ \(t- to)Kh -t 0 ) to <t < 1 1, 


(14) 


(15) 

(16) 

(17) 


and where zo = 2.4 x 10 9 cm, to = 0, and t\ = 3600 s. The 
heating increases linearly in a time much longer than the typical 
dynamic and cooling times (Reale 2010). We set the maximum 
rate Ho = 2 x 10 -3 erg cm -3 s _1 , which can sustain a loop with 
an apex temperature about 3 x 10 6 K according to the scaling 
laws (Rosner et al. 1978). We have checked that the results do 
not change significantly either if we assume a heating function 
that increases exponentially with a comparable rising time, or if 
we consider a magnetic field with half intensity. 

The radial distribution of the heating has a Gaussian shape 
centered on the r = 0 axis. The width of the Gaussian is 
cr H = 3 x 10 8 cm, i.e., 1/10 of the loop’s half-length. The heating 
distribution is uniform with height (z). 


Time [10 3 s] 


M 

0 12 3 



z [1 0 9 cm] 

Fig. 2. Evolution of the density, temperature, and pressure profiles along 
the loop central axis (r = 0) in logarithmic scale after the time- 
dependent heating is switched on. Different line colors mark different 
times, from black (early) to red (late). The profiles are sampled regu- 
larly with a cadence of ^400 s. 


3.1. Results 

The heating is released with gradually increasing intensity 
around the loop’s central axis. Along this axis we expect an evo- 
lution very similar to that obtained from a typical loop model. 
The profiles in Fig. 2 show that this is indeed the case. The tem- 
perature gradually increases in the heated region and reaches 
~3 MK when the heating reaches its maximum rate, i.e., t ~ 
3600 s (red curves in Fig. 2). As is typical of coronal loop mod- 
els, the density increases gradually because of the very slow 
growth rate of the heating, while the plasma evaporates from 
the chromosphere because of the higher pressure driven by the 
heating. The plasma dynamics is minor at all times. The maxi- 
mum apex density (at the final time) is ~2 x 10 9 cm 3 . The loop 
pressure increases uniformly in the corona above 1 dyn/cm 2 . 

Figure 3 shows maps of temperature and density at three dif- 
ferent times, i.e., at the beginning, middle, and end of the heat- 
ing. The heating is conducted only along the magnetic field lines, 
and the proper coronal loop forms and is clearly visible at the 
end of the heating (bottom panel), both in the temperature and 
in the density maps. No local features are visible 1 . Overall, the 
evolution can be described as quasi-static. 

Our attention now turns more specifically to the transition 
region, and in particular to the loop cross section in the expan- 
sion region. Figure 4 shows an enlargement around the region at 
one loop footpoint at two different times. While at early times, 

1 We have verified that, after the heating is switched off, the loop cools 
and drains back to conditions similar to the initial ones. 
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Fig. 3. Density (left, [cm -3 ]) and temperature (right, [K]) in logarithmic 
scale at three subsequent times after the loop heating is switched on: 
beginning (t = 0), half-way (t = 1825 s), and end of the heating (t = 
3610 s). The magnetic field lines are also marked (white lines). 


before the heating is switched on, the transition region layer is 
almost flat with only small ripples, later we see that the transition 
region of the heated part of the flux tube clearly drifts deeper in 
the chromosphere where the heating is more intense, and forms 
a kind of bowl. We now focus on layers at certain fixed temper- 
atures of the transition region, and check whether and to what 


extent the area of that section changes with time because of this 
drift. 

For quantitative analysis, we consider two magnetic field 
lines on the sides of the heated part of the loop. The field lines 
are chosen to bound the width of the heated region in the corona 
(i.e., their distance corresponds to the cr H in the corona). These 
field lines are symmetric and separated from the central axis by 
5 x 10 7 cm at the lower boundary, 10 8 cm in the transition region, 
and 4x 10 8 cm at the loop apex (Fig. 4). They enclose most of the 
heated region (cr H = 3x 10 8 cm in Eq. (15); see Sect. 3). We have 
checked that the results do not change much by choosing other 
field lines on both sides. The field lines experience some read- 
justments during the loop evolution, e.g., they become straighter 
in the low loop region (Fig. 4). This is because the pressure in- 
crease affects the low-beta corona more than the chromosphere 
or transition region, causing it to expand to a greater degree. 

To study the moss visibility we concentrate on layers with 
fixed temperatures, that will be continuously visible in an EUV 
narrow-band detector, for instance. We choose to measure the 
position and width of the layer between the selected field lines 
at temperatures increasing linearly from 10 5 K to 5 x 10 5 K as 
the simulation progresses. Figure 4 shows that both the width 
and the height decrease as the heating increases. As mentioned 
above, the width varies considerably at all selected temperatures 
by -40%, more at higher temperature. This occurs because the 
layers deepen in the chromosphere by a few 10 8 cm, in good 
agreement with the estimate from Eq. (1). The layer becomes 
progressively narrower with time, essentially because the transi- 
tion region drifts deeper inside the chromosphere as the coronal 
pressure increases 2 * in regions where gets smaller. This varia- 
tion is faster at the beginning and then slows down near the end 
of the simulation. The minimum radius of the cross section is 
^0.9 x 10 8 cm and is nearly the same at all selected temperatures 
because the transition region gets thinner as the heating rate in- 
creases. Overall, since the area variation is a tracer of a variation 
in the brightness, this result indicates that even slow changes in 
heating lead to considerable differences in the appearance of the 
moss. 

In both plots at the bottom of Fig. 4, we see periodic oscil- 
lations that modulate the decreasing trend. These are due to the 
heating process. The increasing heating drives moderate pres- 
sure fronts that hit and overshoot in the chromosphere. The pe- 
riod is approximately 300 s which is compatible with the travel- 
ing times of a perturbation at the sound speed (c s - 10 4 a Jf - 
10 7 cm/s) along a tube 5 x 10 9 cm long. 

From the model results, we synthesize the Fe IX line emis- 
sion. In particular, we derive the emission measure in the j-th 
domain cell as em 7 - = n^Vj, where nuj is the hydrogen number 
density in the cell and V) is the cell volume. From the values 
of emission measure and temperature in the cell, we synthesize 
the corresponding emission in the Fe IX line, using CHIANTI. 
Figure 5a shows the 2D map of the emission at time t - 0 in 
the same region as Fig. 4. From the maps at different times, we 
first derive the profiles of the emission along the r-axis at each 
time t by integrating the emission along the z-axis. The time- 
space plot of the emission evolution is then derived from these 
profiles, each normalized to its maximum. The result is shown 
in Fig. 5b. Looking, for instance, at the edge of the red color, we 


2 As long as there are many scale heights in the chromosphere, the 
pressure at the bottom of the chromosphere, Phase, is practically inde- 
pendent of the coronal pressure. The top of the chromosphere (bottom 

of the transition region) occurs at the height z where P = P co r0 na = 

Phase X exp(-z/// g ). 
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t = 0 s 



r [1 0 9 cm] 


t = 3610 s 



-1.0 -0.5 0.0 0.5 1.0 

r [1 0 9 cm] 



Fig. 4. Evolution of the moss region and analysis of the transition region drift and area change: enlargement of density/temperature maps (Fig. 3) 
in the transition region at two times t = 0 and t = 3610 s (see online supplementary movie for an animated version of this evolution). Magnetic 
field lines are marked (white and black lines). The black lines bound the region where we measure the cross section. We also mark the relevant 
cross section with a horizontal segment. The bottom panels show plots of the cross-section radius (left) and of the z position of the layer at five 
linearly increasing temperatures, i.e., 1, 2, 3, 4, 5 x 10 5 K as a function of time, since the beginning of the time-dependent heating. 


can clearly see that the emission shrinks as time progresses and 
the heating increases, and is consistent with the results shown 
in Fig. 4. The oscillations are emphasized because of the depen- 
dence on the square of the density. It can be estimated that the 
radial size is two times smaller. 

4. Discussion 

We investigated an MHD model of an active region loop, and 
in particular the variation in the loop cross section in the moss 
region. To study the low-lying parts of the loops in high detail, 
the assumption of constant cross section all along the loop can- 
not hold, because it is known that the loop expands considerably 
going up from the transition region to the corona (e.g., Gabriel 
1976). We could not help including this effect in the model- 
ing, and how it changes in time and with the plasma (3, and 
this required a proper time-dependent magnetohydrodynamic 
description. 

Our model considers the time-dependent MHD equations in 
a 2D cylindrical geometry including the gravity (for a curved 
loop), the thermal conduction, the coronal heating, and the ra- 
diative losses from optically thin plasma. 

We first created the proper loop topology, with an expanding 
cross section from the chromosphere to the corona, ready for the 


injection of heating. Then, we heated our loop with a large-scale 
slowly-changing (quasi- steady) heating (e.g., Winebarger et al. 
2011; Warren et al. 2010a,b). We focussed on the structure of the 
moss in active regions, i.e., the thin loop layers at temperature 
ranging between about 0.1 MK and 1 MK, that are generally 
interpreted as the footpoints of hotter loops. We measured the 
position and size of the moss layers. 

From Fig. 4 we have seen that the position and size of the 
moss region change with time. Consequently, the emission from 
the moss also varies. 

The primary result from our simulations is that the area of 
the transition region can change substantially if there is a sizable 
change in the quasi-steady heating rate. For example, the area at 
0.5 MK changes by about 40% as the peak temperature in the 
flux tube varies between 1 MK and 4 MK. The plasma beta at 
the relevant heights is low enough that the magnetic field is mini- 
mally influenced by the enhanced gas pressure. The dominant ef- 
fect of the pressure is to force the transition region downward to 
a place where the tube is more highly constricted. This alters the 
DEM(T) curve, since the shape of the curve, including the ratio 
of coronal to transition region emission, depends on how the flux 
tube area varies with temperature. It may not be reasonable to 
independently specify a particular A(T) dependence when mod- 
eling the solar atmosphere. Warren et al. (2010b) found an A(T) 
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Fig. 5. a) Synthetic emission in the Fe IX 171 A line at t = 0, in the 
same enlargement as in Fig. 4. The emission is normalized to its max- 
imum and the color scale is shown in panel b). b) Time-space plot of 
the synthetic emission in the Fe IX 171 A line. We plot the radial pro- 
file of the emission in the vertical direction at all times. Each profile is 
normalized to its maximum. 

that reproduces the observed DEM(T) of an active region core, 
but it is not obvious that this A(T) is consistent with the assumed 
steady heating rate. 

Our model is simplified in several respects. Our 2D cylin- 
drical description assumes a cylindrical symmetry around the 
loop’s central axis. Although there are some indications that this 
might not always be the case (Malanushenko & Schrijver 2013), 
we do not expect considerable differences in the results for mod- 
erate deviations in the shape of the cross section. The structure 
and energy balance of the chromosphere are also simplified into 
an isothermal structure with no radiation transport, but our study 
focusses on layers at much higher temperatures that should not 
be highly affected by this approximation. Finally, although our 
model implies a full MHD description, it is still not entirely self- 
consistent, since the heating is imposed artificially. A fully con- 
sistent model is out of the scope of this work, since we focus on 
the very specific question of the reaction of the low corona to 
heating. 

Another interesting result is the magnitude of the flux tube 
area expansion between the transition region and the corona. The 
area increases by a factor of 2 to 3 over the first 3000 km above 


the chromosphere in our model, depending on the peak temper- 
ature. This increase is much smaller than in the Gabriel (1976), 
Athay (1981), Dowdy et al. (1986) and Warren et al. (2010b) 
models. These models have a throat of rather extreme expan- 
sion, which we suggest may not be realistic. The Gabriel and 
Athay models use a potential field and assume that the flux tube 
is highly constricted at the top of the chromosphere. They also 
assume that the flux tube is isolated, with no surrounding field, 
which is probably a poor assumption, as discussed in Dahlburg 
et al. (2005). Our (3- sensitive MHD model suggests that throats 
of extreme expansion are unlikely, although convective motions, 
not included here, may serve to enhance the effect. We note that 
the dependence of A(T) on heating rate would be even stronger 
if these throats did exist. 
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